More loosely bound hadron molecules at CDF? 
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In a recent paper we have proposed a method to estimate the prompt production cross section 
of A"(3872) at the Tevatron assuming that this particle is a loosely bound molecule of a D and 
a D* meson. Under this hypothesis we find that it is impossible to explain the high prompt 
production cross section found by CDF at cr(X(3872)) ~ 30 ~ 70 nb as our theoretical prediction is 
about 300 times smaller than the measured one. Following our work, Artoisenet and Braaten, have 
£Nj ' suggested that final state interactions in the D°D*° system might be so strong to push the result 

we obtained for the cross section up to the experimental value. Relying on their conclusions we 
show that the production of another very narrow loosely bound molecule, the X s = D a D* s , could be 
similarly enhanced. X s should then be detectable at CDF with a mass of 4080 MeV and a prompt 
production cross section of o~(X a ) ~ 1 -f- 3 nb. 
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Introduction. In a recent paper [1[ we have proposed a method for estimating the prompt production cross 
section of X(3872) [2| at the Tevatron making the assumption that X is a loosely bound molecule of D and D* , 
J^" 1 . with a binding energy as small as So — — 0.25± 0.40 MeV. The motivation for this study is that, after the Belle 
discovery, CDF and DO confirmed the X(3872) in proton-antiproton collisions d, H[ and it seems at odds with 
common intuition that such a loosely bound molecule could be produced promptly (i.e. not from B decay) in a 



high energy hadron collision environment. This was also one of the initial motivations to consider the possibility 
that the X(3872) could be, instead of a molecule, a 'point-like' hadron resulting from the binding of a diquark 
and an antidiquark Q , following the interpretation proposed by Jaffe and Wilczek [§| of pentaquark baryons 
\^ | (antidiquark-antidiquark-quark). 

To start let us summarize the content of [l|. Let us suppose that AT (3872) is an S'-wave bound state of two 
D mesons, namely a l/y/2(D°D* + D°D*°) molecule (we will use the shorthand notation D°D*°) 1 . The 
. molecule production cross section will be proportional to the number of D°D*° pairs in the event. Thus the 
AT (3872) prompt production cross section at the Tevatron could be written as: 
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where k is the relative 3-momentum between the D(pi), D*(y>i) mesons, ipfe) = (X\DD*(k)) is some nor- 
malized bound state wave function characterizing the X(3872). 1Z is the integration region where yj(k) is 
significantly different from zero. The matrix element (DD* (k)\pp) can be computed using standard matrix- 
element /hadronization Monte Carlo programs (MC) like Herwig Q and Pythia [12| 2 . We require our MC tools 
to generate 2 — > 2 QCD events with some loose partonic generation cuts as detailed in [l| 3 . 

As for the determination of the region 1Z in (JIJ we estimate it having in mind a naive gaussian ansatz for 
the bound state wave function. It is straightforward to estimate the momentum spread of the gaussian by 



1 Such a molecule has the correct 1++ quantum numbers of the X(3872). 

2 Open charm meson pairs generated with hadronization Monte Carlo are ordered as a function of their relative center-of-mass 
3-momenta. If more than one D°D*° pair is found in the event, we select the pair having the smaller relative 3-momentum k. 
As a first step we select those pairs which pass the kinematical cuts used in the data analysis made by the CDF collaboration. 

3 Configurations with one gluon recoiling from a cc pair, are those configuration expected to produce two collinear charm quarks 
and in turn collinear open charm mesons. The parton shower algorithms in Herwig and Pythia treat properly these configurations 
at low px whereas they are expected to be less important at higher pj_ . 
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assuming a strong interaction Yukawa potential between the two D mesons. Given that the binding energy £q 
is £ ~ M x -M D - M D , = -0.25 ± 0.40 MeV we find that r ~ 8 fm (8.6 ±1.1 fm) and applying the (minimal) 
uncertainty principle relation, we get the gaussian momentum spread Ap ~ 12 MeV. 

Given the very small binding energy we can estimate k to be as large as k ~ y/2fi(— 0.25 + 0.40) ~ 
17 Me V, n being the r educed mass of D°D*° system, or of the order of the center of mass momentum 
k = y/ \{m Xl m 2 D , m*^) /2mx — 27 MeV. These considerations imply that we can restrict the integration 
region to a ball 71 of radius 4 ~ [0, 35] MeV. 

Keeping k inside 71 we estimate a cr(pp — > AT(3872)) max which is about 30 times smaller than the most 
conservative estimate (a ~ 3.2 nb) of the minimal prompt production cross section measured at CDF [JJ]. In the 
analysis proposed in the experimental cross section is estimated from data to be a ~ 30 4- 70 nb reinforcing 
the negative result obtained with our theoretical calculation that would rather be 300 times smaller than the 
measured one. This fact would undoubtably put in serious trouble the molecular interpretation of V(3872). 

In this note we intend to start from the main result discussed in [l~T| where it is argued that the effect of final 
state interactions in the D°D*° system is such that two corrections should be made to our previous calculation: 
i) the ball 71 should be enlarged to include momenta up to A ~ 300 MeV; ii) a correction factor to the cross 
section we compute (see ([T|)) should be considered so that the actual cross section a* including the full effect of 
final state interaction is 



a*{pp -> V(3872)) = a{k < A) x 
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Assuming that this is the correct way of discussing the V(3872) production we observe that, besides reconciling 
the experimental result with the theoretical computation for the X, this mechanism should enhance the occur- 
rence of an hypothetical new molecule, the D„Dt, which otherwise would be suppressed as one could infer by 
looking at data on D s production at Tevatron [14j (as shown in Fig. [TJ D s is on average ~ 5 times less probable 
than D°). The same data are used to tune our Monte Carlo (Herwig in this calculation) with respect to D s 
production as shown in Fig. [TJ 
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FIG. 1: Differential charm cross section measured in fully hadronic charm decays using 5.8 pb _1 at CDF [TJ. The error 
bars represent the total uncertainty of the measurement. The ratio in the production of D° and D s is 4.4 whereas the 
ratios in the production of D* and D s is about 2.2. The dot-dashed lines are the result of the Monte Carlo simulation 
done with Herwig rescaling the normalizations of the distributions by a factor K = 1.5. This value is in very good 
agreement with the K factor found in [JJ (K — 1.8) using data on da/dA(f>, Acj) being the angle between D° and D*~ 
mesons produced at CDF within some definite cuts in rapidity and transverse momentum. 

The X 8 {1 ++ ) molecule should exist as a partner with strange light quarks of the V(3872). One could expect 
it to be a more compact molecule with respect to the V(3872) as r) particle exchange forces would be at work. 
This enlarges the spread Ap in the relative momentum and naturally makes the ball TZ larger. We will postulate 
a binding energy for X s as small as that of the V(3872) and its mass is expected to be Mx s — 4080 MeV. 
X s could decay into J /ipiriT with a narrow width because of its mass and flavor content (it cannot decay into 
K + K~ J/ip (via <f>) because of phase space; it cannot either decay to J/-0/o(98O) because of quantum numbers), 
or in D s D s j. Using the Herwig hadronization algorithm to compute a(k < A) in (|2"|l we obtain 



a*{pp -> V s (4080)) = 1 4 3 nb 



(3) 



4 Which corresponds to a fco of the Gaussian at ~ 27 MeV and a spread of +12 MeV. 
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FIG. 2: The cross section integrated in bins containing n = 0,1,2,... extra hadrons having a relative momentum 
k < x MeV with respect the D or the D* composing the X(3872) molecule. Following [ll|] we assume that the molecule 
is formed in S-wave with a relative k in the center of mass of D and D* as large as 300 MeV. 



where the value of 3 nb is found pushing the A value up to 600 MeV (following some considerations on the 
possible values of the A cutoff made in [ll|). We obtain definitely similar results using Pythia [l2j . 

Such numbers should put the A" s (4080) molecule in the conditions to be observed at CDF. We would find 
rather surprising that no such state is found assuming that the mechanism ([2]) is correct thus we encourage 
searches of this resonance. 

On the other hand we cast some doubts on the possibility that final state interactions can indeed play such 
a pivotal role as described in [ll[ . First of all we remind that Watson formulae [l3| used in [ll| are valid for 
S'-wave scattering, whereas a relative three-momentum k of 300 MeV indicates that higher partial waves should 
be taken into account. 

Most importantly, we have verified in our MC simulations that as the relative momentum k in the center of 
mass of the molecule is taken to be up to 300 MeV, then other hadrons (on overage more than two) have a 
relative momentum k < 100 MeV with the D or the D* constituting the molecule (see Fig. On the other 
hand the Migdal- Watson theorem for final state interactions requires that only two particles in the final state 
participate to the strong interactions causing them to rescatter. In other words the extra hadrons involved in 
the process do necessarily interfere in an unknown way with the mesons assumed to rescatter into an X(3872). 
This is particularly true as one further enlarges the dimensions of the momentum ball 1Z as required in [11| ■ 

Tetraquarks with a [cs] [cs] might also occur, and one expects the lightest of this family to be a scalar at about 
3930 MeV, as estimated in [15|] . Computing the prompt production cross section is an harder task though. This 
would require some specific model for the fragmentation of partons into diquarks allowing to extract from data 
a ratio of the production rate of [cs] and [cq] diquarks. In turn this would allow, for example, to estimate 
the prompt production cross section of the X s under the hypothesis that the X(3872) produced at CDF is a 
tetraquark. A simple model of parton to diquark fragmentation could be drawn along the lines discussed in (l6j 
where the case of light diquarks was treated. Yet we prefer to postpone such estimate as soon as the first data 
on exotic hadron production will be available from LHCb and ALICE. 

In this note we show that starting from the results discussed in [Tlj we should expect an enhancement in the 
prompt production cross section of an hypothetical new X s (4080) molecular loosely bound resonance constituted 
by a D s D* a pair. We estimate such cross section to be between 1 and 3 nb at the Tevatron. On the other hand 
we cast some doubts on the applicability of the Watson theorem for final state interactions in the calculation at 
hand. We show that in the hadronization shower the number of hadrons in a momentum volume IZ(k) tends to 
grow with k whereas the final state interactions formulae used in [ll[ (see [l3j]) should involve only two hadrons 
at a time. 
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